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Abstract 

In this paper the focus is set on a modified Chua's circuit model equation with saw-tooth function 
in place of piece- wise linear function of Chua's circuit displaying multi-scroll chaotic attractors. 
We study the characteristic properties of first passage times (£fpts) to nth scroll chaotic attractor, 
residence times (£rts) on a scroll attractor and returned times (tRETs) to the middle-scroll attractor. 
tFPTS exhibit a series of Gaussian-like distribution followed by a long tail continuous distribution. 
tRTS and £rets show completely discrete distribution. Power-law variation of mean values of £fpts, 
£rts and £rets with a control parameter is found. On the other hand, mean values of £fpts and 
tRETS have linear dependence with the number of the scroll attractors for fixed values of the control 
parameter. For the system with infinite scroll chaotic attractors normal diffusive motion occurs. 
In the normal diffusion process the mean square displacement grows linearly with time. 

PACS numbers: 02.50.-r, 05.45.-a, 05.60.cd 

Keywords: Modified Chua's circuit equation, multi-scroll attractors, first passage times, residence times, 
returned times, diffusion 
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I. INTRODUCTION 



Chua's circuit is found to show a rich variety of both local and global bifurcations and 
various types of chaotic orbits [Madan, 1993; Tsuneda, 2005; Bilotta et al, 2007]. One of the 
interesting chaotic orbits observed in the Chua's circuit is a double-scroll orbit. Several mod- 
ified Chua's circuit models have been proposed for generating multi-scroll chaotic attractors 
[Lu, & Chen, 2006]. For example, a three-segment piece- wise linear function / characterizing 
volt age- current characteristic of the Chua's diode in the Chua's circuit is replaced by a dis- 
continuous function [Khibnik et al, 1993; Lamarque et al, 1999], sigmoid function [Mahla 
& Badan Palhares, 1993], multiple piece- wise linear segments [Suykens & Vandewalle, 1993; 
Aziz-Alaoui, 1999; Zhong et al, 2002], sign function [Yalcin et al, 2001], sine function [Tang 
et al, 2001], hyperbolic tangent function [Ozoguz et al, 2002], nonlinear term x\x\ [Tang et 
al, 2003], saturated function [Lu et al, 2004] and saw-tooth function [Yu et al, 2007a]. Ex- 
perimental circuit design for generating one-dimensional [Yu et al, 2003], two-dimensional 
[Munoz-Pacheco & Tlelo-Cuautle, 2009] scroll attractors, (2m + 1) x (2m + 1) couple of grid 
multi-scroll attractors [ Yu et al, 2006], hyperchaotic n-scroll attractor [ Yu et al, 2007b] 
and a systematic methodology to design circuits to achieve desired swings, widths, slopes, 
breakpoints, equilibrium points and shapes of n-scroll chaotic attractors [Yu et al, 2005] have 
been proposed. Unidirectionally or diffusively coupling [Kapitaniak & Chua, 1994; Dana et 
al, 2008], cellular neural networks [Arena et al, 1996; Suykens & Chua, 1997], a jerk circuit 
with Josephson junctions [Yalcin, 2007] and thresholding approach [Lu et al, 2008] have 
also been used to generate multi-scroll attractors. In a very recent work Campos-Canton et 
al. [Campos-Canton et al, 2010] constructed a class of three-dimensional dissipative unsta- 
ble systems exhibiting multi-scroll chaotic attractors. Practical applications of multi-scroll 
dynamics are found in broadband signal generators and pseudo-random number generator 
for communication engineering. 

The basic ingredient of multi-scroll attractors in the modified Chua's circuits is the pres- 
ence of multiple breakpoints in the function /. Suppose the number of breakpoints in the 
function / in a multi-scroll Chua's circuit model is very large in a practical consideration and 
theoretically infinite. In such a case f(x) is periodic over the entire range of the argument 
x. The system can exhibit infinite number of scroll chaotic attractors with diffusion charac- 
teristics. In this connection we note that certain nonlinear systems with periodic potential 
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exhibit diffusion dynamics [Blackburn & Gronbech- Jensen, 1996; Popescu et al, 1998; Har- 
ish et al., 2002; Sakthivel & Rajasekar, 2010]. When the mean square displacement, (x 2 ), 
is proportional to t^, /i > the underlying motion is called diffusion, (j, — 1 and /i ^ 1 
correspond to normal and anomalous diffusions. Transition from normal to anomalous and 
vice-versa are found to occur in certain Hamiltonian systems and volume preserving maps 
[Zaslavsky, 2002]. 

Great deal of interest has been focused on the circuit design, identification of different 
forms of the piece-wise linear function f(x) and shapes of various types of scroll attractors 
generated. It is also important to analyse statistical dynamics of multi-scroll chaotic attrac- 
tors and also use nonlinear circuit models to explore statistical properties associated with 
the dynamics exhibited by them. For example, in the multi-scroll chaotic systems, because 
a trajectory escapes from one scroll attractor to nearby scroll attractors, the first passage 
time (FPT), residence time (RT) and returned time (RET) of scroll attractors are impor- 
tant. The focus of the present work is to study the features of these quantities and diffusion 
dynamics in a modified Chua's circuit model equation with multi-scroll chaotic attractors. 
The FPT (£fpt) is the time when a stochastic process x(t) started at t — from a given 
initial value within a domain A of its state space crosses it first time. For systems exhibiting 
multi-scroll chaotic attractors or escape dynamics from one part of an attractor to another 
part of the attractor or from one well of a potential to another well FPT and mean FPT are 
important and study of them has practical applications in many problems [Redner, 2001]. 
RT (£rt) of a scroll attractor is defined as the time duration spend by a trajectory on it 
before passing to another scroll attractor. RET (£ret) of a scroll attractor is the time taken 
by the system's trajectory to re-enter into it. 

We consider the modified Chua's circuit system [Yu et al, 2007] 



x = a[y-f(x)}, 



V = x-y + z, 



z = -py, 



(la) 
(lb) 
(1c) 



where 



f(x) = - £A \ sgn[:r + (2j + I) A] + sgn[x - (2j + I) A] 




3=0 



(Id) 
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with a, p, f , A > 0, N > 1 and 



sgn[xj 



1, if x > 
0, if x = 
-1, if x < 0. 



(le) 



/(x) is a saw-tooth function with amplitude 2A£ and period 2/1. Figure [T] depicts the form 
of f{x) for iV = 2, £ = 0.25 and A = 0.5. The breakpoints a;* of /(x) are indicated by 
painted circles. System (JTJ is capable of generating (2iV + l)-scroll attractors for a range of 
fixed values of the parameters a and /3. In the present work, on the system ([T]) we consider 
the distributions and mean values of FPTs (on the breakpoints x*), residence times (on a 
scroll attractor) and return times (to the middle-scroll attractor) in the case of (2N + 1)- 
scroll chaotic attractors. We study the variation of these quantities with the parameter a for 
N — 1 and with N for few fixed values of a. The probability distribution P(£fpt) exhibits 
a sequence of peaks followed by a long tail continuous distribution. Power-law variation of 
mean FPT with a and linear variation of it with iV for fixed values of a are found, t^rs of 
scroll attractors are discrete. The mean RT also shows power-law dependence on a. For a 
fixed value of a, the mean RT on each attractors are all almost same and independent of N. 
Further, P(£rt) has only a sequence of peaks. On the other hand, P(£ret) of middle-scroll 
attractor has a discrete distribution with number of Gaussian-like envelop. The number 
of such envelops increases with increase in the value of N. The mean RET varies linearly 
with N while it decays following a power-law relation with a. When N = oo (very large 
for a practical consideration) the motion is diffusion along x-direction. The mean square 




FIG. 1: The graph of f(x) given by Eq. JT}i) with N = 2, £ = 0.25 and A = 0.5. The breakpoints 
of f(x) are indicated by painted circles. 
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displacement of the x-component varies linearly with time implying normal diffusion and is 
further confirmed by kurtosis. 



where /' = 

ax 



II. CHARACTERISTICS OF FPT, RT AND RET 

For our numerical study we fix the values of the parameters in Eq. (CQ) as £ = 0.25, A = 0.5, 
(3 = 16 and vary the parameter a. The f(x) curve has discontinuities or breakpoints at 

x* = ±(m-~y m = 1, 2, • • • ,N. (2) 

The equilibrium points of the system ([1]) are given by 

(x e , y e , z e ) = (0, 0, 0), (±m, 0,Tm), m = 1, 2, ■ ■ ■ , N. (3) 

The stability determining eigenvalues of the equilibrium points are the roots of the cubic 
equation 

A 3 + (1 + af')X 2 + (af + (3-a)\ + a(5f = 0, (4) 

= 0.25. For a < a c = 7.74518 all the equilibrium points are stable 
focus with one eigenvalue being real negative and other two being complex conjugate with 
negative real part. As a increases from a small value the magnitude of the negative eigenvalue 
increases while the magnitude of the real part of the complex conjugate eigenvalues decreases. 
At a = a c the two complex conjugate eigenvalues become pure imaginary. For a > a c the 
real part of the these two complex conjugate eigenvalues become positive. So, for a > a c the 
equilibrium points are unstable. In the numerical simulation we found sudden occurrence 
of chaotic motion (crisis) at a = a c . Figure [2] shows the phase portrait of the scroll chaotic 
attractors in x — y plane for a = 10 and for three values of N. 

Now, we present the features of FPT, mean FPT, RT, mean RT, RET and mean RET. 
First we consider the case N = 1. We calculate the FPT on the barrier (breakpoint) x* for a 
set of random initial conditions chosen around the origin for a > a c . Figure |3]^a) shows the 
numerically computed tppT on x* = 0.5 for 5000 initial conditions with a = 10. tppT depends 
on initial condition. £fpts are distributed over a range. Figure El^b) shows the magnification 
of part of Fig. |3]^a). The probability distribution P(tpp^) is calculated using 10 5 £fpts and 
is depicted in Fig. HI P is not a monotonically increasing or decreasing function of FPT. It 
has multiple peaks and a very long tail. The height of the peaks first increases with increase 
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FIG. 2: Phase portrait of (a) three-scroll, (b) five-scroll and (c) eleven-scroll chaotic attractors in 
the x — y plane. The values of the parameters in Eq. ([!]) are a = 10, = 16, £ = 0.25 and A = 0.5. 

in the value of FPT, reaches a maximum value and then slowly decays for larger values of 
FPT. P is continuous for larger values of FPT. For smaller values of £fpt the distribution 
P shows a finite number of Gaussian-like profile. This is clearly seen in Figs.[3](b) and|U(b). 
We note that the peaks of the distribution occur at regular interval of £fpt- Between first 
few bands there exists a range of time interval within which almost no £fpt occurs. The 
above character of £fpt is found for the barrier x* = —0.5 and also for other values of a. 

In Fig. [2] we notice that trajectories leave a scroll attractor through a specific confined 
exit regions. Consider the middle-scroll confined between the breakpoints x* = —0.5 and 
0.5. Through one exit region trajectories enter into left-scroll while through the other they 
enter into right-scroll. A trajectory starting from the vicinity of an exit takes a minimum 
time to arrive at the neighborhood of the next exit. This minimum time is found to be nearly 
unity. Figure [5] shows two trajectories started near the origin and leaving the middle-scroll 
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FIG. 3: (a) Numerically computed FPTs on the barrier x* = 0.5 of the system (pQ) for N = 1 and 
a = 10. (b) Magnification of part of (a) showing finer details. 




FIG. 4: (a) Plot of probability distribution of £fpt on the barrier x* = 0.5 with N = 1 and a = 10. 
(b) Magnification of part of P in (a) corresponding to the Fig. E^b). 

after some time. In this figure for each trajectory there are two vertical lines: the right-line 
denotes the time at which the trajectory leaves the middle-scroll and the left-line marks the 
latest earlier time at which x(t) came near an exit region before leaving the scroll. The time 
difference between the two vertical lines is roughly unity. 

It is interesting to point out the effect of addition of external periodic force on FPT. 
Suppose the system ([1]) is driven by the periodic force gcost. We add this external force 
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FIG. 5: x(t) versus t for two trajectories showing the time of exit from the middle-scroll attractor. 
The breakpoint is x* = 0.5. 

to Eq. (JT^l) . Figure M shows t FPT and -P(£fpt) for two values of g. The continuous and 
nonGaussian distribution observed in Fig. 11(a) for tppT > 12 disappear and a series of 
Gaussian-like distribution occurs even for small values of g. This is evident in Fig. [6(b) 
where g = 0.1. The number of peaks decreases with increase in the value of g and further P 
becomes more and more narrow and the peaks occur at regular interval of time. The time 
intervals between successive peaks of P are almost same and ~ 27r, the period of the external 
force gcost. The periodic nature of P(tppT) can be easily understood. The external force 
periodically oscillates the barrier heights of the function f(x). The chance for crossing the 
critical value x* is large when the barrier height at x* becomes a minimum which happens 
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FIG. 6: 5000 FPTs obtained for the distance x* = 0.5 and the probability distribution of FPT for 
the system ((!]) driven by the periodic force gcost with g = 0.1. 
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once during every period of the external periodic force. This is the reason for the periodic 
occurrence of peaks. If the periodic force is replaced by a Gaussian white noise then the 
FPTs become random, the band-like structure would disappear and P becomes continuous. 
In a purely noise induced process the FPT distribution goes through a single maximum and 
then decays following either exponential or power-law relation. One can expect multi-peaks 
in the presence of both periodic force and weak noise. Suppression of chaotic motion by 
the addition of weak periodic force and delay feedback is realized in several chaotic systems 
[Pyragas, 1992; Braiman & Goldhirsch, 1991]. In the system (JTj) the effect of added periodic 
force is numerically studied for a £ [0, 10] and g £ [0, 0.5]. For a > a c suppression of chaos is 
not observed. For a < a c period-T(= 2ir) orbit is found for the above chosen range of values 
of g. For each fixed value of a < a c , for g value less than a critical value there are three 
coexisting stable periodic orbits. They occur about the equilibrium points of the system. 
As g increases from a small value the size of the orbit increases and above a certain critical 
value the system has a single period-T orbit enclosing all the three equilibrium points. Any 
route to chaotic motion is not found for the parametric choices considered here. 

The mean FPT (£mfpt), the mean time the trajectory initially in the neighbourhood of 
the origin takes to cross the target location x* first time, is calculated for a range of values 
of a > a c . It is obtained averaging over 10 5 FPTs. Numerically computed £mfpt versus a 
is plotted in Fig. [7(a) with x* = 0.5. It decreases rapidly with a following the power-law 
relation: tMFPT = 45.688(a — a c ) -1 ' 003 , i.e., ^mfpt oc l/(ot — a c ). In Fig. EI we can clearly 
see that trajectory moving on a scroll attractor leaves it when enters into a certain interval 
of y called turnstile. We call the length of this interval of y as exit length and denote it 
as y e \. This quantity for the middle-scroll attractor is numerically calculated for a range of 
values of a and the result is shown in Fig. [7J(b). y e \ is found to increase with a — a c obeying 
the power-law relation y e \ = 0.08618(a — a c ) 4127 . The power-law decay of £mfpt with the 
parameter a in Fig. [TJa) can be due to the power-law increase in y e \ with a. 

The characteristics of FPT on x* = 0.5 for N — 1 are found for other values of N and 
x*. In all the cases band-like patterns with gap between the bands for smaller values of 
FPT and continuous distribution for larger values are observed. The range of FPT is found 
to increase with increase in the value of N. This is because the region accessible to the 
trajectories in the region x < increases with increase in the value of N. 

We analyse the dependence of £mfpt (on the target x* = 0.5) with N, the number of the 



10 




FIG. 7: (a) Variation of mean FPT with the parameter a for the target location x* = 0.5. Painted 
circles are numerical data while the continuous line is the best power-law fit. (b) Power-law 
dependence of y e \ with a — a c . 




N 

FIG. 8: MFPT versus N for three fixed values of a. Symbols are numerical data and continuous 
lines are the best straight-line fits. 

saw-tooth function. In Fig. [8] we plotted £mfpt versus N for three fixed values of a. Linear 
variation of ^mfpt with N is observed for each fixed value of a. We obtained 
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14.955iV + 26.46, a = 8.5 



*mfpt(z* = 0.5) = < 12.167iV + 15.753, a = 9 



(5) 



5.996iV + 14.459, a 



10. 



v 



In Fig. |S] we notice that for a fixed value of iV the £mfpt decreases with increase in a. The 
variation of £mfpt with a follows power-law relation while it with N follows linear relation. 
The rate of variation of MFPT with N decreases with increase in a. 

Next, we wish to present the results on residence time statistics. We denote £rt as the 
time the trajectory resides on a scroll attractor before making a transition to any one of 
the adjacent attractors. We fixiV = l and a = 7.52 for which three-scroll chaotic attractor 
occurs. The residence times on the left-, middle- and right-scroll attractors are calculated. 
Unlike £fpts the tjn<s on all the three scroll attractors exhibit a band-like structure over 
the entire possible range of values of t^r- The middle-scroll attractor has more bands than 
the other two scroll attractors. This is because a trajectory can enter into the middle-scroll 
attractor either through the bottom of the attractor (when it enters the attractor from 
the right-scroll attractor) or through the top of the attractor (when it enters the attractor 
from the left-scroll attractor). But a trajectory can enter into the left(right)-scroll attractor 
from the middle-scroll attractor only through the bottom(top) side of the attractor. The 
probability distribution curves shown in Fig. [9] indicates that the bands of RT are very thin. 

When the value of a increases the band-like pattern of RT persists while the maximum 
value of it decreases. Figure ITOl shows mean RT (MRT), £mrt, versus a. Though the number 
of bands of t^x on the middle-scroll attractor is larger than the other scroll attractors the 
MRTs on these three attractors are almost same. For example, for a = 7.47, £mrt for the 
middle and the other two attractors are 393.73 and 389.71 respectively. For a = 7.52 the 
values of £mrt for these attractors are 110.54 and 109.75 respectively. The difference is less 
than 2%. Power- law dependence of MRT on a with the exponent value 1 is observed: 
^mrt = 8.413(a — ct c )~ 0,962 . RTs and MRT on a scroll attractor depend on a but independent 
of N. Thus the MRT on a scroll attractor can be monitored by the control parameter a. 

We believe that the study of RT and MRT on multi-scroll attractors may find practical 
applications. It is noteworthy to mention that RT based detection strategies for nonlinear 
sensors have been proposed [Gammaitoni & Bulsara, 2002; Bulsara et ah, 2003; Dari et 
al., 2010]. This we have explored in the system ([Q). Often spectral technique is used to 
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FIG. 9: Plots of probability distribution of iRT of the three scroll attractors of the system ([I]) for 
N = 1 and ct = 7.52. 




FIG. 10: Variation of tyiRT of the middle-scroll attractor with a for N = 1, tjviRT of the three 
scroll attractors are all almost same. The painted circles are numerical result while the continuous 
line is the best power-law fit. 

detect the presence of weak dc or low-frequency signal. An alternate approach in the case 
of low-frequency signal is the use of tuning of an internal or external noise to induce the 
stochastic resonance. At an optimal noise intensity in a bistable system the signal-to-noise 
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ratio measured at the low-frequency of the signal becomes maximum. To detect a weak 
dc signal one can use residence time asymmetry in noisy bistable devices. The effect of an 
additional dc signal is to skew the potential. In the absence of dc signal the RT distributions 
and MRTs in the two wells of the bistable double-well systems are identical. They will be 
different in the presence of additional dc signal and the difference in the MRT is shown to 
be proportional to the weak dc signal [Gammaitoni & Bulsara, 2002]. This characteristic 
feature can be used to identify the presence of a weak dc signal. In the multi-scroll circuit 
with f(x) given by Eqs. (JTJl-e) the MRTs of a trajectory in all the scrolls are almost identical. 
This property can also be explored for weak dc signal detection as shown below. 
We introduce asymmetry in f(x) given by Eqs. ([Ul-e) by redefining it as 



Figure [TlT a) shows the plot of f(x) versus x for N = 1 with d = and 1. The effect of d 
can be clearly seen. Such forms of f(x) are studied both theoretically and experimentally 
[Aziz-Alaoui, 1999; Lu & Chen, 2006; Suykens & Vandewalle, 1993]. For d = the heights 
of the barriers at the breakpoints x* = 0.5 and —0.5 are same. For d = 1 the barrier height 
at x* = 0.5 remains same while at x* = —0.5 is decreased. This asymmetry in the barrier 
heights can be used to detect the weak signal by measuring the difference in the MRTs on 
the left- and middle-scrolls with and without the dc signal. The MRT of a trajectory in the 
right-scroll remains same for d = and d ^ 0. We define A£mrt = ^mrt(c? = 0) — tuKr(d). 
A^mrt on left- and middle-scrolls are calculated numerically for a = 7.46 and 7.47 for a 
range of values of d. Figure [TlT b) shows the variation of Ai M RT with d. Interestingly, it is 
found to show linear variation with d. From the predetermined relation between A^mrt and 
d we can make an estimate of the value of d by calculating A^mrt- 

We now turn to the discussion on returned time to the middle-scroll attractor (£ret), 
the time taken by a trajectory to re-enter into the middle attractor after leaving it. For a 
trajectory started with an initial condition near the origin 10 5 RETs and the mean value 
of RETs denoted as £mret are numerically calculated for various values of a and N. For 




j=o y - 



where 




d, if x + (2j + I) A < 
0, otherwise. 



(6b) 
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FIG. 11: (a) The graph of f(x) given by Eq. © with N = 1, £ = 0.25, A = 0.5 and for d = and 
1. (b) Dependence of A^mrt = ^mrt(^ = 0) — £mrt(^) of left - and middle-scroll chaotic attractors 
on d for two values of a. Curves 1 and 2 are for middle-scroll attractor with a = 7.46 and 7.47 
respectively. Curves 3 and 4 are for left-scroll attractor with a = 7.46 and 7.47 respectively. The 
symbols are numerical data and continuous lines are the best straight-line fit. 

N = 1 three discrete bands of RETs occur. The number of bands increases with increase 
in the value of N. P^ret) obtained with 10 5 RETs for N = 1, 2 and 3 are reported in 
Fig. [T2l Distribution of RETs is different from those of FPT and RT. The width of all the 
bands of RETs are very small. In the Figs. IT27 b) and \V2lc) we can clearly notice a series 
of Gaussian-like profile of P. The first profile is predominant over the others. The height 
of the profiles decays rapidly with £ret- The above characteristics of RETs are found for 
higher values of N also. 

Next, we observe different kinds of dependence of ^mret on iV and a. In Figs. fTBTa) we 
plotted £mret versus iV for various fixed values of a. For each fixed value of a the £mret 
increases linearly with N. This is because the number of scroll attractors increases linearly 
with N. We obtained the following linear fits: 

a = 7.8 : £ M ret = 47.920iV - 28.602. 
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FIG. 12: Probability distribution of returned times to the middle attractor for three values of iV 
with a = 9. 




FIG. 13: Variation of mean returned time to the middle attractor with (a) N and (b) a. In (a) the 
values of a for the top to bottom curves are 7.8, 8, 8.5, 9 and 10 respectively. In (b) the values of 
N for the bottom to top curves are 1, 2, 3, 4 and 5 respectively. The symbols are numerical data 
while the continuous lines are the best fits. 



16 



a = 8.0 : t M RET = 30.243iV - 18.995. 

a = 8.5 : £ M ret = 18.557iV - 9.381. 

a = 9.0 : t M RET = 12.639iV - 7.926. 

a = 10.0 : £ M ret = 7.329iV - 1.668. 
The rate of divergence of £mret with N decreases with increase in a. Regarding the depen- 
dence of tMRET on a we notice power-law variation in Fig. ITBT b). The following result is 
obtained: 



N 


= 1 : 


^MRET 


= 9.499(a - 


- a c )~ 


-0.923 


N 


= 2 : 


^MRET 


= 27.055(a 


- a c ) 


,-0.913 


N 


= 3 : 


^MRET 


= 45.618(a 


~ «c) 


,—0.922 


N 


= 4 : 


^MRET 


= 67.093(a 




,-0.980 


N 


= 5 : 


^MRET 


= 85.992(a 


- a c ) 


,-0.993 



III. NORMAL DIFFUSION 

In this section we report the occurrence of normal diffusion in Eq. ([1]). When N = oo the 
function f(x) has infinite number of saw-tooth segments and for a > a c the motion is not 
bounded to a finite range of x in the limit t — > oo. Figure IT41 shows x versus t for a = 8 and 
9. For a = 8 the divergence of x is very slow. We can notice relatively a higher divergence 
of x for a = 9 in Figs. [T4T c) and IT4T d). The motion is essentially diffusive. 

In order to capture the type of diffusion we calculate the mean square displacement over 
a set of M initial conditions given by 

M 

(x\t))= hm ^E^W-^W)] 2 , (7) 

M— s-oo i = \ 
t— ¥oo 

where x^ 1 ' is the ith. trajectory and (x^(t)) is the mean value of x^'(t). In our numerical 
simulation M = 10 4 and the initial conditions are chosen around the origin. The variation 
of (x 2 (t)) against t is shown in Fig. [15] for three values of a. On the log 10 — log 10 scale the 
slopes of the best straight-line fits for a = 8, 9 and 10 are ~ 1. In other words, (x 2 (t)) ~ 
as t — > oo with ji — 1. 

A linear time dependence of (x 2 (t)) is the hallmark of normal diffusion and such a diffu- 
sion process can be described by a classical random walk. The value of the exponent /i is 
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FIG. 14: x(t) versus t for two values of a of the system (JTJ) when N = oo. The subplots (b) and (d) 
are magnification of detail evolution of x(t) of part of the solution shown in (a) and (c) respectively. 

calculated for a range of values of a in the interval [a c , 11]. In all the cases \i is found to be 
1 and anomalous diffusion (/i ^ 1) is not found. We note that anomalous diffusion is found 
in several conservative systems when the phase space has accelerating modes or stochastic 
layers. In these systems when a trajectory comes closer to such an accelerating mode or a 
stochastic layer the system is accelerated rapidly and (x 2 ) diverges nonlinearly with time 
during the trapped times and moves chaotically otherwise. The combination of these two 
types of motion leads to anomalous diffusion. In certain dissipative chaotic diffusive systems 
divergence of x takes place either during laminar intervals or chaotic evolution of a state 
variable. In the system ([T]) there is no laminar phase. The growth of the state variable x is 
due to the passage of the trajectory to the adjacent scroll attractors along one direction and 



18 




FIG. 15: Plot of (x 2 (t)} as a function of time in log 10 — log 10 scale for (a) a = 8, (b) a = 9 and 
(c) a = 10 for the system ([1]). The symbols are numerical data and continuous lines are the best 
straight-line fits. 

also because of the presence of infinite number of scroll attractors. This is clearly evident 
in Fig. fl4T d). In this way diffusion in the system ([1]) is different from the diffusion in the 
systems like damped and forced pendulum. 

The coefficient D c in (x 2 (t)) = D c t increases with increase in the value of a. The values 
of D c for a = 8, 9 and 10 are 0.018, 0.146 and 0.346 respectively. For a one-dimensional 
piece- wise linear periodic map Schuster and Just [Schuster & Just, 2005] shown that the 
presence of diffusion implies chaotic motion and further obtained the scaling law associated 
with the diffusion coefficient D c . For the conservative standard map an analytical expression 
for momentum distribution at time n for momentum scales larger than the control parameter 
K is obtained [Ott, 1993]. Damped oscillatory variation of D c with the control parameter 
K about the quasilinear value of D c has been reported [Rechester & White, 1980]. 

Normal diffusion can be characterized by kurtosis which is defined as 

{{x-{x)Y) 

((x-(x)) 2 ) 2 ' K) 

For normal diffusion K = 3 in the limit t — > oo. In Fig. [16] the variation of K is plotted for 
three values of a. For large t the value of K is ~ 3 confirming normal diffusion. 

We studied the occurrence of diffusion in two other multi-scroll chaotic systems, one with 



f(x) = J2^L-n(~ l) J_1 tanh[fc(x — 2j)\ [Ozoguz et al., 2002; Salama et al., 2003], hyperbolic 
tangent function, realizable using a multi-cycle transconductors composed of alternating D- 
P cells and another one proposed by Tang and his co-workers [Tang et al., 2001] in which 
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FIG. 16: Kurtosis K versus t for a = 8 (marked by painted circles), 9 (open circles) and 10 (painted 
triangles) for the system (pQ). 

f(x) is a sine function. Normal diffusion is found in both the systems. 

IV. CONCLUSION 

In the present work we have considered the modified Chua's circuit model equation 
[Eq. (JTJ)] with saw-tooth function for f(x) proposed to generate multi-scroll chaotic at- 
tractors. Attention has been paid to the study of characteristics of FPTs, RTs and RTs. 
We have presented the distribution of these quantities and explored the dependence of mean 
values of them on the control parameter a and the number of scroll attractors. From our 
study we can make some general results for infinitely large number of scroll attractors. Be- 
cause the forms of the function f(x) associated with all the scroll attractors are identical 
the resident time characteristics of all the scroll attractors are almost the same within the 
numerical accuracy. Further, the returned times statistics of all the scroll attractors are also 
same. In the system with the infinite number of scroll attractors the motion is shown to be 
normal diffusion type. We wish to point out that the multi-scroll systems and the features of 
multi-scroll attractors can be useful to emulate various logic gates and the ability to switch 
easily between the different operational roles. In a recent work Murlai et al. [Murali et al., 
2009] have shown that the response of a bistable system with a two square waves and appro- 
priate noise intensity produces a logic outputs (NOR/OR) and (NAND/AND). We believe 
that similar results can be realized if one replaces the bistable system by the multi-scroll 
systems. Further the different scrolls can be used to represent different logic states in which 
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case MRT on a scroll attractor is important. In the system ([I]) MRT on a scroll can be 
varied over a wide range of the control parameter a as shown in Fig. [TUl 
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